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It is a firm prediction of the concordance Cold Dark Matter (CDM) cosmological model that 
galaxy clusters live at the intersection of large-scale structure filaments^ The thread-like 
structure of this "cosmic web" has been traced by galaxy redshift surveys for decades^^. 
More recently the Warm-Hot Intergalactic Medium (WHIM) residing in low redshift fila- 
ments has been observed in emission'* and absorptionPll However, a refiable direct detection 
of the underlying Dark Matter skeleton, which should contain more than half of all mattei^, 
remained elusive, as earUer candidates for such detections''^^" were either falsifiediiii^ or 
suffered from low signal-to-noise ratios* and unphysical misalignements of dark and lumi- 
nous matter^ Here we report the detection of a dark matter filament connecting the two 
main components of the Abell 222/223 supercluster system from its weak gravitational lens- 
ing signal, both in a non-parametric mass reconstruction and in parametric model fits. This 
filament is coincident with an overdensity of galaxies^" and diffuse, soft X-ray emission^ 
and contributes mass comparable to that of an additional galaxy cluster to the total mass of 
the supercluster. Combined with X-ray observations'*, we place an upper limit of 0.09 on the 
hot gas fraction, the mass of X-ray emitting gas divided by the total mass, in the filament. 

Abell 222 and Abell 223, the latter a double galaxy cluster in itself, form a supercluster 
system of three galaxy clusters at a redshift of 2 ~ 0.21 [ref. T3^, separated on the sky by ~ 14'. 
Gravitational lensing distorts the images of faint background galaxies as their light passes massive 
foreground structures. The foreground mass and its distribution can be deduced from measuring 
the shear field imprinted on the shapes of the background galaxies. Additional information on this 
process is given in the supplementary information. The mass reconstruction in Figure [T] shows 
a mass bridge connecting A 222 and the southern component of A 223 (A 223-S) at the 4. la 
significance level. This mass reconstruction does not assume any model or physical prior on the 
mass distribution. 



1 



To show that the mass bridge extending between A 222 and A 223 is not caused by the 
overlap of the cluster halos but in fact due to additional mass, we also fit parametric models to the 
three clusters plus a filament component. The clusters were modelled as elliptical Navarro-Frenk- 
White (NFW) profiles'"* with a fixed mass-concentration relation'^. We used a simple model for 
the filament, with a flat ridge line connecting the clusters, exponential cut-offs at the filament end 
points in the clusters, and a King profile'-^ describing the radial density distribution, as suggested 
by previous studie^i^III. We show in the supplementary information that the exact ellipticity has 
little impact on the significance of the filament. 

The best fit parameters of this model were determined with a Monte-Carlo Markov Chain 
(MCMC) and are shown in Fig.[2j The likelihood-ratio test prefers models with a filament compo- 
nent with 96.0% confidence over a fit with three NFW halos only. A small degeneracy exists in the 
model between the strength of the filament and the virial radii of A 222 and A 223-S. The fitting 
procedure tries to keep the total amount of mass in the supercluster system constant at the level 
indicated by the observed reduced shear. Thus, it is not necessarily the case that sample points 
with a positive filament contribution indeed have more mass in the filament area than a 3 clusters 
only model has. The reason is that the additional filament mass might be compensated for with 
lower cluster masses. We find that the integrated surface mass density along the filament ridge line 
exceeds that of the clusters only model in 98.5% of all sample points. This indicates that the data 
strongly prefers models with additional mass between A 222 and A 223-S and that this preference 
is stronger than the confidence level derived from the likelihood-ratio test. The difference is prob- 
ably due to the oversimplified model, which is not a good representation of the true filament shape. 
The data on the other hand is not able to constrain more complex models. Extensions to the simple 
model we tried were replacing the flat ridge line with a parabola as well as replacing the King 
profile with a cored profile leaving the exponent free. The latter was essentially unconstrained. 
The parabolic ridge line model produced a marginally better fit that, however, was statistically 
consistent with the flat model. Moreover, the likelihood-ratio test did not find a preference for the 
parabolic shape. 

The virial masses inferred from the MCMC are lower than those reported earlier for this 
system'", which were obtained from fitting a circular two-component NFW model to A 222 and 
A 223. Compared to this approach, our more complex model removes mass from the individual 
supercluster constituents and redistributes it to the filament component. Reproducing the two- 
component fit with free concentration parameters, which was used in the previous study, we find 
M2oo(A 222) = (2.7lo:7) x 10^^ Mq, which is in good agreement, and M2oo(A 223) = (3.4l}:|j) x 
10''^ Mq, which overlaps the Icr error bars of the earlier study. Here and in the following, all error 
bars are single standard deviations. 

The detection of a filament with a dimensionless surface mass density of k ~ 0.03 is un- 
expected. Simulations generally predict the surface mass density of filaments to be much lowei'^ 
and not to be detectable individually'**. These predictions, however, are based on the assumption 
that the longer axis of the filament is aligned with the plane of the sky and that we look through 
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the filament along its minor axis. If the filament were inclined with respect to the line-of- sight and 
we were to look almost along its major axis, the projected mass could reach the observed level. A 
timing argument'^ ^" can be made to show that the latter scenario is more plausible in the A 222/3 
system. In this argument we treat A 223 as a single cluster and neglect the filament component, 
such that we have to deal only with two bodies, A 222 and A 223. The redshifts of A 222 and 
A 223 differ by A2; = 0.005, corresponding to a line-of-sight separation of 18 Mpc if the redshift 
difference is entirely due to Hubble flow. Let us assume for a moment that the difference is caused 
only by peculiar velocities. Then at 2 = 00, the clusters were at the same location in the Hubble 
flow. We let them move away from each other with some velocity and inclination angle with re- 
spect to the line-of-sight and later turn around and approach each other. The parameter space of 
total system mass and inclination angle that reproduces the observed configuration at z = 0.21 is 
completely degenerate. Nevertheless, in order to explain the observed configuration purely with 
peculiar velocity, this model requires a minimum mass of (2.61 ± 0.05) x 10^^ Mq with an incli- 
nation angle of 46 degrees, where the error on the mass is caused solely by the uncertainty of the 
Hubble constant. Since this is more than 10 standard deviations above our mass estimate for the 
sum of both clusters, we infer that at least part of the observed redshift difference is due to Hubble 
flow, and that we are looking along the filament's major axis. 

The combination of our weak-lensing detection with the observed X-ray emission of 0.91 ± 
0.25 keV WHIM plasma'* lets us constrain the hot gas fraction in the filament. Assuming that the 
distribution of the hot plasma is uniform and adopting a metallicity of Z = 0.2 Solar, the mass 
of the X-ray emitting gas inside a cylindrical region with radius 330 kpc centred on (01:37:45.00, 
12:54:19.6, Figure |3]) with a length along our line-of-sight of / = 18 Mpc, as suggested by our 
timing argument, is Mgas = 5.8 x 10^^ Mq. The assumption of uniform density is certainly a strong 
simplification. Because the X-ray emissivity depends on the average of the squared gas density, 
a non-uniform density distribution can lead to strong changes in the X-ray luminosity. Thus, if 
the filament consists of denser clumps embedded into lower density gas (as has been observed in 
the outskirts of the Perseus Cluster^^), or even if there is a smooth non-negligible density gradient 
within the region used for spectral extraction, then our best fit mean density will be overestimated. 
The quoted gas mass should therefore be considered as an upper limit, and the true mass may be 
lower by up to a factor of 2-3. 

We estimated the total mass of the filament from the reconstructed surface mass-density map 
and the model fits within the same region where we measured the gas mass. The conversion of 
dimensionless surface mass density to physical units requires knowledge of the source redshifts. 
We randomly sampled galaxies with our Rc-band magnitude distribution from photometric redshift 
catalogues'^. The mean redshift of these random catalogues is Zs = 1.2. We emphasize that for 
a cluster at ^ = 0.21, the error in mass caused by the uncertainty of the redshift distribution 
is small. An error as large as A^;^ = 0.2 causes only a 5% error. In the reconstructed K-map, 
the mass inside the extraction circle is Mgi = (6.5 ±0.1) x 10^'^ Mq, where the error is small 
due to the highly correlated noise of the smoothed shear field inside the extraction aperture. For 
the parametric model fit, the inferred mass is higher but consistent within one standard deviation. 
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Mfil = (9.8±4.4) X 10^^ Mq. The corresponding upper limits on the hot gas fractions vary between 
fx = 0.06 — 0.09, a value that is lower than the gas fraction in galaxy clusters^^. This is consistent 
with the expectation that a significant fraction of the WHIM in filaments is too cold to emit X-rays 
detectable by XMM-Newton^'*. 
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Figure 1 : Mass reconstruction of A 222/223. The background image is a three colour- 
composite SuprimeCam image based on observations with the 8 m Subaru telescope 
during the nights of Oct. 15 (A 222) and 20 (A 223), 2001 in V-, Rc- and i'-bands. 
We obtained the data from the SIVIOKA science archive (http: //smoka.nao . ac. jp/). 
The FWHIVI of the stellar point-spread function varies between QH57 and (/.'TO in our final 
co-added images. Overlayed are the reconstructed surface mass density (blue) above 
K = 0.0077, corresponding to S = 2.36 x 10^^ IVIpc^, and significance contours above 
the mean of the field edge, rising in steps of 0.5a and starting from 2.5a. Dashed con- 
tours mark underdense regions at the same significance levels. Supplementary Figure 1 
shows the corresponding B-mode map. The reconstruction is based on 40,341 galax- 
ies whose colours are not consistent with early type galaxies at the cluster redshift. The 
shear field was smoothed with a 2' Gaussian. The significance was assessed from the 
variance of 800 mass maps created from catalogues with randomised background galaxy 
orientation. We measured the shapes of these galaxies primarily in the Rc-band, sup- 
plementing the galaxy shape catalogue with measurements from the other two bands for 
galaxies for which no shapes could be measured in the Rc-band, to estimate the gravita- 
tional sheari5lll6]. a 222 is detected at ~ 8.0a in the south, A 223 is the double-peaked 
structure in the north seen at ~ 7a. 
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Figure 2: Posterior probability distributions for cluster virial radii and fila- 
ment strength. Shown are the 68% and 95% confidence intervals on the cluster 
virial radii 

'"2oo( ) SiUd the filament strength kq. The confidence intervals are derived 
from 30,000 IVICIVIC sample points. The filament model is described by K{6,r) = 

Kq |i + exp [{\6\ - Oi)/a] + (r/rc)^| , where the coordinate 6 runs along the filament ridge 
line and r is orthogonal to it. This model predicts the surface mass density at discrete grid 
points from which we computed our observable, the reduced shear, via a convolution in 
Fourier space. The data cannot constrain the steepness of the exponential cut-off at the 
filament endpoints a and the radial core scale r^. These were fixed at their approximate 
best-fit values of a = 0.45 IVIpc and Vc = 0.54 IVIpc. The data also cannot constrain the 
cluster ellipticity and orientation. These were held fixed at the values measured from the 
isodensity contours of early-type galaxies^^. The ratios of minor/major axes and the po- 
sition angles of the ellipses are (0.63, 0.69, 0.70) and (65°, 34°, 3°) for A 222, A 223-S, and 
A 223-N, respectively. We further explore the impact of cluster ellipticity on the filament 
detection in the supplementary information. 
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Figure 3: Surface mass density of the best fit parametric model. The surface mass 
density distribution of the best parameters in Fig. [2] was smoothed with a 2' Gaussian 
to have the same physical resolution as the mass reconstruction in Fig. [T| The yellow 
crosses marl< the end points of the filament model. These were determined from the 
visual impression of the filament axis in Fig.[l] The IVICIVIC is not able to constrain their 
location. In the model, the filament ridge line is not aligned with the axis connecting the 
centers of A 222 and A 223-S. This is a fairly common occurrence (~ 9%) for straight 
filament but may also indicate some curvature, which occurs in ~ 53% of all intercluster 
filaments^'' and is not included in our simple model. Overlayed are X-ray contours from 
XIVIIVI-Newton observations'^ (red) and significance contours of the colour-selected early- 
type galaxy density^" (beige), showing the alignment of all three filament constituents. 
The blacl^ circle marl<s the region inside which the gas mass and the filament mass were 
estimated. 



Supplementary Information 

In this document we provide additional information on the data reduction, shear measurement, and 
tests performed. 

SI Data 

We reduced the data using the publicly available theli pipelin^. The data were astrometri- 
cally calibrated against the USNO-Bl.O catalogue, which has a precision of ~ 0'.'2. Since this is 
insufficient to align objects on dithered exposure without introducing errors in the gravitational 
shear measurement, theli uses the SWarp calibration software to align images internally. The 
astrometric residuals in this process are smaller than 0'.'02 or a tenth of a pixel. 

No photometric standards were observed in the nights our data was taken. We calibrated our 
Rc-band data against the catalogue of our earlier study of this system'^. We fixed the zero-points 
of the other two bands by matching the expected colours of stars from the Pickles stellar librar)^^ 
to the observed colours of stars in the field. 

The gravitational shear was estimated independently in all three passbands using the lensfit 
algorithml2lIll!. Lensfit is a Bayesian model fitting code, which fits the sum of an exponential and 
a de Vaucouleurs profile to galaxies. The model fit has six free parameters, the galaxy position 
{x,y), its ellipticity (61,62), brightness, bulge-fraction, and galaxy scale r. The model fitting is 
done in Fourier space so that the brightness and position can be marginalised over analytically. 
The radius and bulge-fraction are marginalised over numerically, leaving a likelihood surface as a 
function of (61, 62). Lensfit acts on the individual exposures and generates a posterior probability 
for each galaxy by summing over the posterior ellipticity distributions generated by measurements 
in each exposure. 

Because lensfit is a forward fitting method, its model space needs to be convolved with a PSF 
model describing the stellar point spread function at the location of each galaxy. We modelled the 
spatial variation of the PSF across the focal plane with a bivariate polynomial, whose zeroth and 
first order coefficients were allowed to vary from chip to chip. 

Galaxy shapes are highly correlated across different passbands and simply combining shear 
estimates from different passbands entails the risk of degrading them to the band with the worst 
seeing or the largest remaining systematics. To maximise the lensing information from all three 
passbands we therefore chose a different route. We used shape measurements from the Rc-band 
data, which has the best seeing and is the deepest. For galaxies without Rc-band shape information, 
we used i'-band measurements, and in turn V-band measurements were used for galaxies that also 
have no i'-band shapes. 



We used an empirical approach to find an optimal balance between the rejection of clus- 
ter/foreground galaxies and the density of our lensing catalogue. We placed aperture mass filterP 
on the three known structures in the field (A 222, A 223-S, A 223-N). For galaxies with two 
colours measured, we rejected galaxies inside a parallelogram placed on a colour-colour diagram. 
An MCMC algorithm moved the vertexes of this parallelogram to maximise the aperture mass 
signal. A similar procedure was used for galaxies with one colour measured, either V — Rc or 
Rc — i'. A rectangle inside which galaxies were rejected from the lensing catalogue was placed 
on a colour-magnitude diagram. Again, the corners of this rectangle were moved to maximise the 
aperture mass signal. For both procedures, a maximum magnitude, down to which the colour cut 
was applied, was determined at the same time. Finally, for galaxies observed in only one passband, 
we imposed a magnitude cut in the same way. 

This method rejected galaxies with colours consistent with early-type galaxies at the cluster 
redshift and showed that fainter galaxies produce a stronger aperture mass signal than brighter 
ones. Both properties are indications that we measure a real cluster signal and that the redshift 
scaling of the lensing signal is at least approximately correct. At the same time we find that colour 
and magnitude cuts made to maximise the aperture mass signal of the known clusters, also lead to 
an increase in the lensing strength of the filament. This supports our view that the filament is a real 
structure located at the same redshift as the galaxy clusters. 

We also note that the filament is seen in mass reconstructions generated from lensing cata- 
logues of the individual passbands, but at lower levels consistent with the shallower data. In the 
Rc-band reconstruction, a mass bridge is present at the 3.5a level. In the V- and i'-band it is seen at 
the 2a and 2.5a level, respectively. We emphasise again that one cannot simply combine these sig- 
nificances in quadrature because a large covariance exists between the shear estimates in different 
passbands. 

Finally, we show E- and B-mode reconstructions of our combined lensing catalogue to test 



for PSF modelling systematics (Figure S 1 1 



S2 Fitting elliptical profiles 

The ellipticity of the clusters connected by a filament has the potential to be confused with the 
filament itself and may in extreme cases even lead to the visual impression of a filament even if 
the clusters' mass distribution is completely described by elliptical NFW halos. In this section we 
further explore the influence of cluster ellipticity on the significance of our filament detection. 

To recapitulate, we modelled the clusters with fixed ellipticities measured from the distribu- 
tion of their early type galaxies. In this configuration, A 223-S points almost exactly along the 
filament ridgeline and A 222 is about halfway turned towards the filament. Consequentially both 
clusters contribute more mass towards the filament area than they would in a spherical configu- 



ration. This is reflected in small differences in the filament significance if we fit three spherical 
clusters. In this case, the integrated mass along the filament exceeds that of the three clusters only 
case in 98.9% compared to 98.5% for the fixed ellipticities. The likelihood-ratio test prefers three 
spherical clusters plus a filament over a model with only three spherical clusters at 97.8% confi- 
dence, compared to the 96.0% significance for the fixed ellipticities used in the main paper. The 
filament mass in the case of spherical halos is higher by about 30%. This is within the error bars 
we report but easily explained by the diminished contribution of A 222 and A 223-S to the surface 
mass density in the filament region. 

We used fixed cluster ellipticities because the data is not able to constrain the ellipticity, as 
evidenced by Figure |S2j However, even though we allowed for more elongated halo shapes and 
for A 222 to point exactly along the filament ridgeline, the marginalised probability distribution 
of Kq shows that a positive filament distribution is still strongly favoured. Again the integrated 
surface mass density along the filament exceeds that of a model without filament component for 
98.4% of all sample points and the likelihood ratio test prefers the presence of a filament at 98.2% 
significance. The latter number is noticeably higher than what we report for fixed cluster elliptici- 
ties. The reason is that the "best fit" values, which are really unconstrained, have A 223-S oriented 
perpendicular to the filament. As a result, the presence of a filament is even more strongly required 
than in the configuration we considered in the main body of the paper. The likelihood-ratio test 
shows that elliptical clusters are slightly preferred by the data over spherical ones (84.3%), but - 
as suggested by the flat posteriors on g - no preference exists for the specific values of axis ratios 
and position angles chosen by us (19.5%). In light of what the data tell us, we thus consider our 
choice of cluster ellipticities to be conservative. 
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Figure S1 : E- and B-mode reconstructions. Panel a again shows the mass recon- 
struction from Figure [1] The colour scale in the background is dimensionless surface 
mass density and the contours are significance contours starting at 2(j and rising in steps 
of la above the mean of the field edges (solid b\acW), excluding the top edge, which is 
strongly affected by stray-light. Dashed white contours are at the same negative levels. 
Panel b shows a mass reconstruction of the B-modes, obtained by rotating all galaxies 
by 45°. IVIore and/or higher peal^s than expected from a Gaussian random field would 
be an indicator of systematic residuals in the shear estimation. The peak count here is 
consistent with expectations for a pure noise field^°. The shear field in both panels was 
smoothed with a 2' Gaussian. Both panels have identical colour scales. 
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Figure S2: Posterior probability distributions with free ellipticities. Similar to Fig- 
ure [2] we show the 68% and 95% confidence contours on the cluster virial radii r2oo(-) 
and filament strength kq. In addition we tried to constrain the cluster axis ratios g( ) and 
position angles 0( ) for A 222 and A 223-S. Their values for A 223-N were held fixed at 
those used in the main paper, as the mass distribution of A 223-N shows no degeneracy 
with the filament strength. We put a flat prior 0.5 < g < 1 on the axis ratio to encompass 
all but the most elliptical halos^^ This also includes the values used in the main paper. 
Allowing for more extreme axis ratios does not seem to be warranted by the absence of a 
trend in the filament strength with axis ratios. 



